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Abstract. 

We explore long-range interactions between two atoms excited into high principal 
quantum number n Rydberg states, and present calculated potential energy curves for 
various symmetries of doubly excited ns and np rubidium atoms. We show that the 
potential curves for these symmetries exhibit deep GHz) potential wells, which can 
support very extended (~ //m) bound vibrational states {macrodimers). We present n- 
scaling relations for both the depth of the wells and the equilibrium separations Re 
of these macrodimers, and explore their response to small electric fields and stability 
with respect to predissociation. Finally, we present a scheme to form and study these 
macrodimers via photoassociation, and show how one can probe the various ^-character 
of the potential wells. 



PACS numbers: 32.80.Rm, 03.67.Lx, 32.80.Pj, 34.20.Cf 
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1. Introduction 

Rydberg atoms have long been studied because of tlieir peculiar properties such 
as long lifetimes, large cross sections, and very large polarizabilities These 
exaggerated properties lead to strong interactions between the Rydberg atoms, which 
have been experimentally detected in recent years {2!, ^ ■ Such strong Rydberg- Rydberg 
interactions have fueled a growing interest in the field of quantum computing, and 
over the past decade, their application for quantum information processing, such as 
fast quantum gates jj, [sf, or quantum random walks ^ have been proposed. Also 
of particular interest is the excitation blockade effect j7|, where one Rydberg atom 
actually prevents the excitation of other nearby atoms in an ultracold sample 



This phenomenon was recently observed in microtraps [isl, and a C-NOT gate was 



implemented using the behavior 15 



Another active area of research with Rydberg atoms is the predicted existence of 
long-range "exotic molecules". In one scenario, one atom remains in its ground state, 
while another atom is excited to a Rydberg state. The most famous examples of this type 
of interaction are the trilobite and butterfly states, so-called because of the resemblence 
of their respective wave functions to these creatures. The theoretical framework for such 



interactions was first proposed in [16|, but were not observed until more recently in [1 



The second type of long-range interaction is predicted to occur when both atoms are 



excited to Rydberg atoms. In 18|, it was first predicted that weakly bound macrodimers 



could be formed from the induced Van der Waals interactions of two such excited atoms. 



However, more recent work [19[ has shown that larger, more stable macrodimers can be 
formed from the strong mixing between ^-characters of various Rydberg states. Recent 
measurements have shown signatures of such macrodimers in spectra of cesium Rydberg 



samples [20 



In this article, we present long-range potential energy curves corresponding to the 
interaction between pairs of rubidium atoms excited to ns and np Rydberg states. 
In general, Rydberg- Rydberg interactions will only mix states that share the same 



molecular symmetry [21|, |22|. Thus, only common symmetries between the excited 
Rydberg molecular state and the state to which it is most strongly coupled are 
relevant. For rubidium, the doubly excited ns atom pair is most strongly coupled 
to the np+ [n — l)p asymptote, while the doubly excited np atom pair is most strongly 
coupled to the ns + {n + l)s asymptote. Since all ss' states have rrij = ±|, the only 
common symmetries with any pp' state are Q = |mj^ + mj^l = 0, 1. In this manner, we 
find that the relevant symmetries for the doubly excited ns and np asymptotes are 0^, 
0~ and lu- We analyze all three cases for both pairs and show that potential wells exist 
for all of them. We also describe in detail properties of the bound levels within each 
well. 

The paper is arranged as follows: in §2l we review how to build the basis states 
used to compute the potential energy curves at long-range, and describe the existence of 
potential wells for certain asymptotes. In ^|3l we investigate the effects of small external 
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electric fields on the potential curves, and in §U we discuss the scaling of the wells with 
principal quantum number n. Finally, in ^ we calculate bound levels supported in 
those wells, and estimate their lifetimes. We also outline how photoassociation could be 
used to form and probe macrodimers. This is followed by concluding remarks in ^ 



2. Molecular Curves 

2.1. Basis States 

In this section, we review the general theory for calculating the interaction potential 
curves. These curves are calculated by diagonalizing the interaction Hamiltonian in the 
Hund's case (c) basis set, which is appropriate when the spin-orbit coupling becomes 
significant and fine structure cannot be ignored, as is the case here. 

We first consider two free Rydberg atoms in states \a) = \n,i,j,mj) and 
\a') = \n',i',j',m'j), where n is the principal quantum number, £ the orbital angular 
momentum, and rrij is the projection of the total angular momentum j = i + s onto 
a quantization axis (chosen in the 2;-direction for convenience). The long-range Hund's 
case (c) basis states are constructed as follows: 

|a; a'; Qg/u) ~ |a)i|a')2 - p(-l)^+^'|a')i|a)2, (1) 

where Q = rrij + m'j is the projection of the total angular momentum on the molecular 
axis and is conserved. The quantum number p describes the symmetry property under 
inversion and is 1(— 1) for g[u) states. 

For Q = 0, we need to additionally account for the reflection through a plane 
containing the internuclear axis. Such a reflection will either leave the wave function 
unaffected or it will change the sign of the wave function. We distinguish between these 
symmetric and antisymmetric states under the reflection operator as follows: 

K/u)=^-^\%/u), (2) 
where behaves according to the following rules 



23 



24 



a^\A) =(-1)^1 -A) (3) 
a,\S,Ms) = {-lf~^''\S,-Ms). (4) 

References 22| and 25| give the technical details for determining which states 
comprise the basis of the np + np rubidium asymptote. Although the procedure to find 
the basis states for different molecular asymptotes, such as ns+ns, np+np, or nd+nd, is 
the same, the states making up these basis sets, in general, will be different. We do not 
review the procedure for building the basis sets here, but we note that all relevant {i.e. 
strongly coupled) molecular asymptotes within the vicinity of the asymptotic doubly- 
excited Rydberg state being considered are included in each respective basis set. We 
again note here that because doubly excited ns (np) rubidium atoms are most strongly 
coupled to pp' {ss') states and because only common symmetries of such Rydberg states 
are allowed to mix, the relevant symmetries for the ns + ns and np + np asymptotic 
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Table 1. Asymptotic 0+ molecular states included in the Rb 70p + 70p basis set, 
which diagonalize the interaction Hamiltonian (see text). The basis states have been 
symmetrized with respect to the reflection operator ([2]) and each \ai; a2',0g) state is 
defined by equation ([T]). 
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Rydberg states that we consider are 0^, 0„ and 1„. As an example, Table [T] lists the 
basis set for the 0^ symmetry near the Rb 70p + 70p molecular asymptote. 

2.2. Long-range Interactions 

The interaction matrix we consider consists of both the long-range Rydberg-Rydberg 
interaction and the atomic fine structure. Here, "long-range" refers to the case where 
no electron exchange takes place i.e. the electronic clouds about both nuclei do not 
overlap. This occurs when the distance R between the two nuclei is greater than the 



LeRoy Radius [26 



Rlr = 2 [{n,h\r'\n,hy/' + {n.i^lr'ln.i,)'/'] . (5) 

When the distance between the two atoms is larger than Rlr, the interaction between 
them is described by the residual Coulomb potential between two non-overlapping 



charge distributions |27|, which can be truncated to give only the dipole-dipole (V^) 
and quadrupole-quadrupole (V^) terms. For two atoms lying along the 2;-axis, the Vd 
and Vq terms can be simplified to give (28| : 

Vl{R) = -tD^ltt 2^ B^rYrih)Y,-{h). (6) 
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is the binomial 



Here, L = 1(2) for dipolar (quadrupolar) interactions, 5^ = j^,^ "^^^| 
coefficient, rj is the position of electron i from its center, and L = 2L + 1. 

Because the molecular basis states are linear combinations of the atomic states 
determined through symmetry considerations, each matrix element will actually be a 
sum of multiple interactions, i. e. 

{a;a';Qg/u\VL\b;b';Qg/^) = {a; a'\VL\b;b') - Pai-lY^+'H^'; ci\VL\b;b') 

-p,{-lY>'+'i{a;a'\VL\b';b) 
+ pM-lY^+''^+'>'+'^{a';a\VL\b';b) , (7) 

where \a;a') = \a)i\a')2 and so on. An analytical expression for the long-range 
interactions is obtained using angular momentum algebra in terms of 3-j and 6-j 
symbols: 

-T-.L -T-.L 
'^13 '^24 



(8) 




where jtot = Y.t=iji^ = + = 2ji + 1, 
element. For |1; 2) 



and TZfj = {i\r^\j) is the radial matrix 
1 3; 4), the matrix element is given by (j8]) plus the sum of the two 
atomic asymptotic Rydberg energy values. That is: 

{l;2\VLiR)\l;2) = m + Ei + E2, (9) 

where rii is the principal quantum number and Si is 
). Since Ai = dipole transitions are 



29 



and 



30 



with Ei given by - — — 

2{ni - diY 

the quantum defect (values given in 
forbidden, only the L = 2 term of equation ([8]) will contribute in ([9]). 

Figure [1] shows the results of diagonalization for the 0^, 0~ and 1„ symmetries of 
the 70s + 70s, and 70p + 70p asymptotes with no background electric field. In all of 
these plots, the energies are measured from the ionization threshold of rubidium. We 
see that all three symmetries feature large potential wells and for the remainder of this 
paper, we focus our attention on formation properties of bound states within these wells 
and analyze the stability of these macrodimers. We note here that we also explored the 
potential curves near the nd + nd asymptotes, but no wells were found; hence we do not 
display those curves here. 



3. Electric Field Dependence 

Production and/or detection of macrodimers will rely on external electric fields. In 
addition, since experiments cannot completely shield the atoms from undesired stray 
fields, it is important to study their effect on our calculated curves. 
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Figure 1. (Color online) Long-range interaction curves for the 0,„ (left), 0+ (middle), 
and 1„ (right) symmetries of doubly-excited ns (top row), and np (bottom row) 
Rb Rydberg atoms near n = 70. We highlight the potential wells and label their 
corresponding asymptotic state. Note: All energies are measured from the Rb 
ionization threshold. 



Strictly speaking, applying an external electric field F breaks the D^oh symmetry of 
homonuclear dimers, and consequently, the basis states defined by (II]) would no longer 
be valid. In principle, one then needs to diagonalize the interaction matrix in a basis 
set containing every possible Stark state, as was done in 3l| . However, since the effects 
of such an electric field should be adiabatic, we assume that the Dooh symmetry is still 
approximately valid for small electric fields. 

We consider the effects of such an electric field as a perturbation to the original 
Hamiltonian. In general, an applied electric field will define a quantization axis; 
the molecular axes of our macrodimers will then be at some random angle to this 
quantization axis. In that case, one needs to transform the molecular-fixed frame back 
into the laboratory- fixed frame 32| . To simplify our calculations, we assume that the two 
Rydberg atoms are first confined in an optical lattice, such that the quantization axes 
of the macrodimer and the electric field coincide (see figure [2]). Such one-dimensional 
optical lattices have already been used to experimentally excite Rydberg atoms from 



small Bose-Einstein condensates located at individual sites 33 . We envision a similar 
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Figure 2. (Color online) Two Rydberg atoms confined within harmonic traps provided 
by optical lattice sites. The electric field F is directed along the z-axis, coinciding with 
the molecular axis of the dimer. We also indicate in this cartoon that the two excitation 
lasers (see section r5.2[) propagate in the y-direction, but the polarization directions of 
both lasers are along the z-axis. 



one- dimensional optical lattice with the distance between adjacent (or subsequent) sites 
adjusted to coincide with Re, but containing a single atom per site. The lattice could 
be switched off during the Rydberg excitation to allow a cleaner signal. 

For an electric field directed along the ^-axis, the perturbation Hamiltonian for a 
single atom is given by Fr cos 9, where F is the magnitude of the electric field, r is the 
distance of the valence electron from its nucleus, and 6 is the angle between F and r. 
We express the new eigenstates (Stark or dressed states) of the perturbed Hamiltonian 
as a series expansion using the unperturbed asymptotic eigenstates, i.e. 

\a)k = J2^k,^{F)\a,)k, (10) 

i 

where \ai)k = \ni,ii,ji,m^j^)k are the unperturbed (undressed) atomic states of atom k, 



and bk,i{F) are field-dependent eigenvectors, resulting from diagonalization [34|. Here, 
the index i stands for nj, ii, and jj. In Fig. |3l we show the Stark map for \me\ =0 near 
n = 69. Although the limits of this summation are technically rii — Umax (where rimax 
is the highest n value in the basis) and ii — (n^ — 1), we restrict the summation to 
{n — 2) < rii < {n + 2) and ii <3] the bi{F) coefficients are insignificant (> 2 orders of 
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Figure 3. Atomic Rydberg energies vs. the electric field for Rb with |m^| = near 
n = 69: the curves labeled by n = 69 include states with i> 3. 



magnitude smaller) for states lying outside these bounds. Using the dressed states (ITOll 
in ([T]), we define the dressed molecular states as: 

\a)i\a)2 = 5^6i,.(F)62,i(F)|a..)i|a,)2 • (11) 

We then use this basis to redefine the properly symmetrized dressed molecular basis 
given in Table [1] and to diagonalize the Rydberg- Rydberg interaction matrix. 

In Fig. m we illustrate the effect of F on the curves near 70p + 70p of the 0^ 
symmetry, in a side-by-side comparison of the curves for F = and 0.3 V/cm. The 



atomic Stark states were computed using the method described in [18|, |3J] and the 
interaction curves for the "pseudo-symmetry" were obtained by diagonalization of the 
Rydberg- Rydberg interaction matrix in the dressed molecular basis set f lTT]) . We find 
some minor differences: the Stark effect is most notable in the shifting of the potential 
curves, especially the asymptotic energies. However, the relative shape of the curves are 
only slighty changed and most importantly, the large potential well is robust against 
small electric fields. 



4. Scaling 

We focus our attention on the wells correlated to the (n — l)s + (n + 2)s asymptote 
near np3/2 + np3/2, and those correlated to (n — l)pi/2 + ^^^1/2 near ns + ns, for F = 0. 
We calculated the curves for a large range of n, and found that the wells follow simple 
n-scaling behaviors (see Fig. |5]). As will be discussed in ^ the wells are produced by 
the mixing of several states with different ^-character. The exact mixing takes place 
mainly via dipole and quadrupole interactions, and occurs between nearby n-states: 
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Separation R (1000 a^) Separation R (1000 a^) 

Figure 4. (Color online) 0+ molecular curves for Rb 70p + 70p: (a) F = 0, and (b) 
0.3 V/cm. The states involved in the well correlated to 69s + 72s are identified (see 
section [5]). The zcro-cncrgy is set at the 70p3/2 + 70p3/2 asymptote in both plots. 



hence the wells depend on the combination of multipole interaction between states, and 
the proximity (in energy) of those states. 

To derive a simple ra-scaling behavior for both the depth Dg and equilibrium 
separation Re, we assume that the dipole-dipole coupling is the dominant interaction 
between states, and that the wells are formed as a result of an avoided crossing between 
two potential curves (see Fig. |5]^a)). As mentioned above, the real situation is much 
more complex, but these assumptions allow for a simple treatment. 

The energy difference AEa/s = Ea — Ef^ is defined by the difference between the 
asymptotes of the two crossing states, a and /3. Here, a = Ua^s + rta^s-, with energy 



2E„ 



and (3 = rip-^p + rii^^p, with energy 2Ei3 



In both 



energy terms, z/^. = n^. — 5^. includes the principal quantum number n and quantum 



defect 5 of the appropriate atomic state 7i of atom i. Assuming that the relevant atomic 
states in a given asymptote are separated by An of the order unity, we can expand the 
energies as 



so that 



2E^ 



AE^ 



and 2Efi 



(12) 



+ 



Anp 



v. 



AUr 



v. 



Arir. 



Av 



where we assume v^^ ~ z/q,^ + Av with Av of order one, and A = {Ana — Arip] 



(13) 
and 



B = SAvAng. In the cases leading to our wells. Aria — Au/^ so that the leading 
dependence is AEa/3 oc n~^. 




Figure 5. (Color online) (a) Isolated avoided crossing in the 0~ symmetry curves 
of doubly excited ns atoms, which results in the potential well correlated to the 
69pi/2 + 70pi/2 asymptote. We assume the interactions at the crossing arc mostly 
dipolar in nature (see text) and label the well depth D,, and the equilibrium separation 
Re at the avoided crossing. We also demonstrate that at long range the curves behave 
as ~ i/R^, but in the i?-range of the avoided crossing, they behave as ~ 1/R^. (b) 
Scaling relations for the well depth De vs n for the 0~ symmetry of ns + ns (top) and 
np + np (bottom), (c) Scaling relations for the equilibrium separation R,, vs n for the 
same curves as in (b). 



/^From the sketch depicted in Fig. [5]^a), assuming leading dipole-dipole interactions, 
the equihbrium separation Re occurs at the "intersection" of two attractive and 
repulsive curves separated by AE^is, i.e., AEajs — Ca/Rl — +Ci3/Iil, which lead to 
i?e ~ [{Cp + Ca\/ AE^pY^'^ . Our assumption that the two crossing curves behave as 
~ is valid in the region of the intersection; at larger values of i?, however {R > 

80 000 Oq), the curves behave more like ~ 1/R^ (see figure [5]). /^From the scaling 
C/3 + C*o oc and AE^jB oc ra"^, we obtain R^ oc n^/'^ . As for the dissociation energy 
-De, it is simply given by AEajS — Ca/Rl — ^e, which scales as AE^p, i.e. oc n~'^. 
Fig. El^b) shows a plot of vs. n for the 0~ symmetry of the ns + ns and np + np 
asymptotes, indicating that Dg indeed scales more like ~ n~^ + (blue curve) than 
purely n~^ (red curve). For the same wells, Fig.[5](c) shows that Re follows the predicted 
n^^'^ scaling. In the interest of space, we do not show plots for all of the symmetries 
highlighted in Fig. [H but we note that we find the same approximate scaling for all 
other symmetries. 
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Although the analytical derivations above give good agreements with numerically 
determined values of Dg and Re, the slight discrepencies, especially with the np + np 
plots, reflect the more complex nature of the interactions. For example, quadrupole 
coupling is present in our calculations (although its effect is generally small). We also 
point out that in the three np + np cases, the formation of each well is not clearly given 
by an avoided crossing of two curves, but rather by several interacting curves (see next 
section). Nonetheless, the good agreement depicted in Fig. [5] indicate that these more 
complicated interactions act only as small corrections. 

5. Forming Macrodimers 
5.1. Energy Levels and Lifetimes 

The wells identifled in Fig. [1] support many bound levels. We list the lowest levels for 
each well in Table [21 together with the corresponding classical inner and outer turning 
points. For those wells around the ~ 70 asymptotes, the deepest energy levels are 
separated by about 1—2 MHz, corresponding to oscillation periods of a few fis, rapid 
enough to allow for several oscillations during the lifetime of the Rydberg atoms (roughly 
a few hundred /is for n = 70). These energy splittings also allow for detection through 
spectroscopic means. As illustrated by the values of the turning points in Table [21 the 
bound levels are very extended, leading to macrodimers of a few fim in size. 

As described in ^|21 the molecular curves are a direct result of the ^-mixing 
occuring between the electronic basis states ([1]): each molecular electronic state \xx{R)) 
(corresponding to the potential curve Ux{R)) is expanded onto the electronic basis states, 
the amount of mixing varying with R: 

\xx{r)) = J24'\r)\j)^ (14) 

where c^j''^\r) are the eigenvectors after diagonalization for each separation R, and |j) 
the electronic basis states ([T]). 

In [loj , we showed that the potential well corresponding to the 69s + 72s asymptote 
of the 0^ symmetry curves near Rb np + np was composed of flve nearby asymptotes. 
We do not review the detailed treatment here, but highlight the signiflcant asymptotes 
in Fig. M^a) (left panel). We note that the composition of this well is due to the strong 
dipole mixing between the flve highlighted states, which lead to the well having not 
only ss' character, but also pp' character. In the right panel of Fig. M^a), we illustrate 
the same information for the well correlated to the 69pi/2 + 70pi/2 asymptote near the 
70s + 70s asymptote of the 0^ symmetry. Again, we flnd that the states contributing 
the most to the existence of this well correspond to asymptotes above 69pi/2 + 70pi/2. 

In Fig. Mjo), we depict the ^-mixing leading to the potential wells in Fig. M^a): 
this is given by the \cj{R)\'^ coefficients. For both wells near np + np and ns + ns, the 
molecular basis states |j) mixed by the dipole and quadrupole interactions correspond 
to asymptotes that lie above the asymptote of the wells. In the case of 70p + 70p, 
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Table 2. Energies of the six deepest bound levels (measured from the bottom of the 
well) and corresponding classical turning points for the O^J", 0~ and symmetries near 



doubly excited ns and np Rb Rydberg atoms near n ~ 70. 


Asymptote Symmetry 


V 


Energy (MHz) 


Ri (a.u.) 


i?2 (a.u.) 


ns + ns 0+ 
y 





1.035 


46,137 


46,538 




1 


3.121 


45,985 


46,688 




2 


5.353 


45,870 


46,800 




3 


7.477 


45,780 


46,888 




4 


9.567 


45,702 


46,963 




5 


11.645 


45,630 


47,032 


ns + ns 0^ 





1.034 


45,072 


45,488 




1 


3.161 


44,913 


45,650 




2 


5.262 


44,803 


45,761 




3 


7.331 


44,716 


45,854 




4 


9.368 


44,639 


45,934 
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Figure 6. (Color online) (a) 0+ symmetry curves of the rubidium 70p + 70p molecular 
asymptote (left panel), and the 70s + 70s molecular asymptote; both plots are zeroed 
at the ionization level of rubidium. For both panels, we highlight the molecular curves 
corresponding to the five electronic states contributing the most to the formation 
of the well (see text), (b) Composition of the 69s+72s well (left panel) and of the 
69pi/2 + 70pi/2 well (right panel): probabilities |cj(i?)p of the electronic states that 
contribute the most to the formation of the well vs. the nuclear distance R. Inset: 
zoom of the inner region. 



the 69s + 72s molecular level couples strongly to both the 69p + 71p states (above) 
and the 68p + 72p states (below). However, the relative energy differences between the 
asymptotes results in a much stronger interaction between 69s + 72s and the 69p + 71j9 
states than with the 68p + 72p states. This is why there is little contribution from 
the 68p + 72p states in the formation of the well. In the case of ns + ns, the states 
directly below the well correlated to the 69pi/2 + 71pi/2 asymptote are 68p + 71p states. 
In general, the strength of quadrupole coupling between np states and n'p states is 
very weak. Combining this with the large spacing between the asymptotic energy levels 
results in minimal contributions from the 68p + 71p states in the formation of the well. 

The insignificant contributions of the states below the potential wells in both 
cases indicate little chance of predissociation to these lower asymptotes, and hence the 
macrodimers should be long-lived (limited only by the lifetime of the Rydberg atoms). 



In [19[, we showed that the nonadiabatic coupling between the 69s + 72s curve and 
the curves immediately below were very small, leading to metastable macrodimers. For 
the wells near the ns + ns asymptotes, we reach the same conclusion, i. e. metastable 
macrodimers with lifetime limited by that of the Rydberg atoms. 



5.2. Photoassociation 

Exciting two ground state atoms into a bound level via photoassociation (PA) will 
allow us to probe the different electronic characters mentioned in 15.11 In the following 
treatment, we describe a PA scheme for the formation of macrodimers bound by the 
highlighted well of the 0^ symmetry of the np + np asymptote (see Fig. [6|^a) left panel). 
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We assume that the ground state atoms are first excited to intermediate Rydberg 
states (treated as the "ground" states) so that the coupling to higher Rydberg states 
is enhanced. Since the bound states of this particular well have electronic character 
that is mostly \ns;n's) and \np;n'p), we assumed two possible "ground" states making 
transitions to different ^-character in the well: 40p3/2 + 40p3/2 to the ss' components, 
and 41s + 41s to the pp' components. For simplicity, we choose intermediate states near 
n ~ 40 (see inset of figure [7]) because the electronic potential curves of these states are 



asymptotically flat in the R region of the potential well we wish to populate (see [19[); 
this greatly simplifies the calculations for the PA rates. However, our calculations can 
of course be modified to fit other experimental parameters and conditions. 

Figure [7] shows a schematic for our proposed two-photon formation mechanism for 
the case of doubly-excited np atoms. The macrodimers we predict could be realized 
and identified spectroscopically by red-detuning the excitation lasers from the resonant 
frequency of the 70p3/2 -l- 70p3/2 molecular Rydberg level. 



The PA rate K„ for two atoms into a bound level v can be calculated [35| using 

K„ oc /1/2 |(0^,|(XA|eVlr2|x5)|0^,)|^ (15) 

where Ii and I2 are the intensities of laser 1 and 2, |0„(i?)) and \xx{R)) are the radial 
and electronic wave functions inside the well, respectively, \(j)g{R)) and \xg{R)) are the 
radial and electronic wave functions of the ground state, respectively, and rj and e are 
the location and charge of the electron i. Using the expression (IHj) for \xx{R)), and 
assuming that \xg) is independent of R (corresponding to a fiat curve), we can rewrite 
( IT5|) as 

/■oo 2 

(16) 



poo 

K„cx JiJ2y;i(rfit/2),f / dR<p:iR)c*iR)ijg{R) 

J ^0 



where di = {njij\eri\ngig) and ^2 = {n'ji'j\er2\n'gi'g) are the electronic dipole moments 
between electronic states ja^; a'g) and \aj; a'j) for atom 1 and atom 2, respectively. 

In 19|, we presented PA rate calculations for bound levels in the potential well 
correlated to the 69s -|- 72s asymptote from a flat radial ground state distribution. The 
results of the PA rate against the detuning A from the atomic 70p3/2 levels are shown 
in Fig. [Ht^a). Since the general expression for the PA rate f lT6|) is proportional to the 
dipole moments and the laser intensities, our calculated rates are given in arbitrary 
units set to a maximum of one (for the strongest rate starting from 41s + 41s). In fact, 
once a pair of atoms has been excited to the "ground" state with fixed laser intensities, 
the PA rate plotted in Fig. |8] represents the probability of forming a macrodimer; the 
transition from the 5s Rb atoms to the intermediate "ground" state atoms can easily 
be saturated so that the PA process always starts with a pair of 40p3/2 + 40p3/2 or 
41s + 41s. We compare these PA rates to those obtained if the ground state radial wave 
function (j)g{R) is assumed to be a gaussian centered on R^. with a standard deviation 
of 14 500 ao (roughly half the FWHM of the potential well) in Fig. [H^b). The choice of 
a gaussian approximates the wave function obtained by thermally averaging harmonic 
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Figure 7. (Color online) Our proposed two-photon pliotoassociation scheme for the 
formation of 70p + 70p rubidium Rydberg macrodimers. The ground state atoms are 
populated to a bound level inside the well by dual lasers, each of which is red-detuned 
from the resonance signal of the 70p3/2 -I- 70p3/2 molecular Rydberg state. Inset: Each 
atom is initially excited to an intermediate Rydberg state, which is considered to be 
the "ground" state in our discussion (see text). We note that the single atom detuning 
levels are not to scale. 



oscillator wave functions over a harmonic trapping potential [e.g. in an optical lattice) 
for both "ground" state atoms. In those plots, the PA rates starting from both atoms 
in 40p3/2 are shown in red and from 41s in turquoise, respectively. The rapid oscillation 
between large rates for an even bound level {v = 0,2,4, . . .) and small rates for odd 
levels {v = 1,3,5, . . .) gives the apparent envelope of the PA signal. This behavior is 
due to the oscillatory nature of the radial wave functions inside the well </)„(-R): the 
integral in equation ( 1T6|) will be near zero for odd wavefunctions. We highlight this for 
the gaussian ipgi^R) distribution of 40^3/2 + 40p3/2 in (c) and 41s -|-41s in (d), where we 
show a zoom of the deepest levels v (on a log-scale). 

We see in plots (a) and (b) of Fig. [8] that the signature of a macrodimer would 
manifest itself by the appearance of a signal starting at A ~ —0.93 GHz red-detuned 
from the 70p3/2 atomic level, and ending abruptly at ~ —1.36 GHz. The shape of the 
signals indicate that the rates can reveal details of the ^-mixing inside the potential 
well; for example, both plots show that the K„ mimic the probabilities \cj{R)\'^ shown in 
figureint^b). The progressive decrease of the 41s signal beginning at —1.36 GHz, followed 
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Figure 8. (Color online) PA rate vs. the detuning A from the atomic state 

for (a) a constant radial ground state distribution and (b) a gaussian radial ground 
state distribution - both on linear scales set to a maximum of 1. The |41s;41s) state 
populates the pp' character in the well (turquoise), and the |40p3/2; 40p3/2) the ss' 
(red). Plots (c) and (d) show the rates for the deepest levels of the gaussian distribution 
on a logarithmic scale (see text). 



by sharp increases between —1.22 and —1.16 GHz correspond to the slow decreases of 
the pp' components between i? ~ 40 000 — 50 000 ao and their sharp increases around 
-R ~ 33 000 — 35 000 ao- For the 40p3/2 signal, the major feature common to both 
(a) and (b) is the significant drop in K„ between —1.13 and —1.09 GHz, which mirrors 
the decrease in the ss' states between i? ~ 52 000 — 55 000 in Fig. [6]^b). As noted 



m 



19[, this range of frequency with a noticeable drop in the PA rate could serve as 
a switch to excite or not excite a macrodimer, depending on the "ground" state being 
used. We also note that in both cases, the signals for the 40p3/2 + 40p3/2 "ground" state 
are higher overall across the R regime of the well, with a few exceptions. This indicates 
that despite the presence of the pp' character, this well is still largely composed of ss' 
character. 

As expected, the gaussian ground state distribution shares many qualitative 
features with the constant ground state distribution. We considered both ground state 
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distributions having populations in the R range of ~ 30 000 to 70 000 ao- Normahzing 
both ground state distributions over this range yielded slightly larger signal rates in 
the deepest part of the well for the gaussian, corresponding to its peak. However, the 
major difference between Fig. [8]^a) and (b) is noticeable in both the 40p3/2 + 40^3/2 and 
the 41s + 41s rate signals between A ~ —1.09 and —0.85 GHz. Whereas the uniform 
distribution shows a steady 41s + 41s signal and a steady increase of the 40p3/2 + 40^3/2 
signal, the gaussian distribution shows both signals rapidly decreasing. As we chose to 
center the gaussian at the minimum of the well, the decrease in both signals obviously 
corresponds to the decreasing probability of the gaussian distribution in this R regime 
[i.e. the tails). If one wanted to take advantage of the large, isolated 69s72s character 
at higher R, this could easily be accomplished by recentering the gaussian appropriately. 



6. Conclusion 

We have presented long-range interaction potential curves for the 0^, 0~, and 1„ 
symmetries of doubly excited ns and np Rydberg atoms and have demonstrated the 
existence of potential wells between these excited atoms. These wells are very deep 
and very extended, due to the strong ^-mixing between the various electronic Rydberg 
states. These wells are robust against small electric fields and support several bound 
vibrational levels, separated by a few MHz, which could be detected in spectroscopy 
experiments. The macrodimers corresponding to these bound vibrational levels are 
stable with respect to predissociation and have lifetimes limited only by the Rydberg 
atoms themselves. These macrodimers could be realized through population of the 
vibrational energy levels by photoassociation, resulting in a detectable signal that could 
be used to probe the various ^-character of the potential well. 

In conclusion, we note that the detection of such extended dimers could facilitate 
studies in a variety of areas. For example, the effect of retardation on the interaction 
at very large separation, which becomes important if the photon time-of-flight between 
the atoms is comparable to the classical orbital period of a Rydberg electron around 
its core ISj, could potentially be probed experimentally. Another example relates 
to chemistry of molecules with high internal energy; a third atom approaching a 
macrodimer could quench its internal state to lower levels, or could potentially react 
with the molecule at very large distances and create a new product such as trilobite or 



butterfly Rydberg molecules [16|. Finally, as mentioned in the introduction, Rydberg 
atoms are being investigated intensively for quantum information processing, e.g. using 
the blockade mechanism [3], and the possibility of frequency ranges where the PA 
rate is strong or weak due to the ^-character mixing could potentially be used as a 
quantum mechanical switch. These few examples illustrate some possible applications 
of macrodimers. 
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